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TITLE 

MEASURING METHOD FOR DECIDING DIRECTION TO A FLICKERING 
SOURCE 

5 

TECHNICAL FIELD 

The present invention relates to a method for deciding 
the direction to a flickering source in relation to a 
measurement point in an electrical network with 
10 alternating current with a network frequency with low- 
frequency amplitude variations from a flickering source. 
The present invention also relates to an arrangement 
| comprising mcano for carrying out the method. 

15 BACKGROUND ART 

During the production of electricity, generators are used 
that produce an alternating voltage around a certain 
frequency. The user of the produced electricity is aware 
of which frequency is involved and which voltage is 

20 provided on the electricity network. The user of the 
electricity network requires as clean electricity as 
| possible^ — that — That is , the user requires electricity 
that is well-defined at the specified frequency and the 
specified voltage. In normal electricity networks that 

25 supply towns and factories, it is, however, normal for 
the electricity network to comprise low- frequency voltage 
distortions. Voltage distortions arise when various 
devices are cyclically connected into the network, which 
devices comprise, for example, capacitances and 

30 inductances. 



Large power-hungry installations such as induction 
| furnaces, compressors, lift motors, pumps, etc_;_, help act 
to increase the amount of voltage distortion in the 
electricity network. These loads often cause low- 
5 frequency (<25 Hz) and periodic fluctuations in the 
effective value of the voltage. The phenomenon is usually 
called flicker. The most noticeable effect of flicker is 
the annoyance that can be experienced when the light 
intensity from incandescent lamps fluctuates in line with 

10 the voltage variations caused by the flicker. Studies 
have shown that people are particularly sensitive to 
light fluctuations with a repetition frequency in the 
range 0.5 Hz to 25 Hz. At the level of maximum 
sensitivity (approx. 9 Hz), the relative voltage change 

15 needs only to be 0.25% for people to be able to perceive 
that the light from the incandescent lamp flickers. 
Problems with flicker arise primarily in areas with very 
j heavy industry (ironworks and paper-mills, etc^) but can 
also arise in areas with weak electricity networks and in 

20 the vicinity of wind-power installations. 



In the event of a dirty electricity network comprising 
flicker, it is of interest to know where the flickering 
source is located. An electricity producer wants to show 
25 that it is a consumer that makes the electricity network 
dirty and can, in such a case, demand that the consumer 
pay a penalty charge or puts — te-h-e — problem — right corrects 
the problem to eliminate the flicker . A consumer, on the 
other hand, wants of course to show that it is not him 
the consumer causing the problem. In addition, the 
consumer wants to have value for money and thus wants to 
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show that it is the electricity producer that is 
supplying dirty electricity. 



j Putting — right Correcting the problem of flicker is often 
5 an expensive operation where all or parts of the network 
can need to be reconstructed in such a way that the 
internal impedance of the electricity network is reduced, 
for example by moans — use of new and larger cables. 
Another way to — correct the problem right — is to 

10 install counter measures against the flicker in question. 
Such counter measures are usually very expensive. An 
example of a counter measure is "Static Var Compensation 
(SVC)" which dynamically controls changes in the system. 

15 Methods exist for determining the content of voltage 
flicker, which are described in the standard IEC-61000-4- 
15. Measuring instruments that record voltage flicker 
according to this standard show the presence of voltage 
flicker by calculating and displaying the parameters Ifl, 

20 | Pst and Pit. There i.Sj_ however^ nothing in the standard 
or in existing measuring instruments that shows the 
direction to the flickering source in relation to the 
measurement point. 

25 There is thus a great need for a method and an 
arrangement that can determine the direction of 
propagation of flicker. The direction of propagation of 
the flicker shows whether the flickering source is above 
or below a measurement point and can thus be of great use 

30 when the flickering source is to be traced. 

DISCLOSURE OF INVENTION 
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The present invention intends to solve the problems that 
are described above of low-frequency flicker that gives 
rise to periodic fluctuations in the effective value of 
the voltage. In the following, such periodic fluctuations 
5 are called low-frequency variations or flicker. 



the direction to a flickering source in relation to a 
measurement point in an electricity network with 

10 alternating current with a network frequency f c with low- 
frequency amplitude variations from the flickering 
| source. The method is — characterized — ift provides that an 
amplitude-modulated current signal and an amplitude- 
modulated voltage signal are recorded. Both the 

15 amplitude-modulated current signal and the amplitude- 
modulated voltage signal are signal processed in such a 
way that the low-frequency amplitude variations in both 
the current and voltage are separated from the carrier 
wave in the form of a flicker component for current and a 

20 flicker component for voltage. The flicker component for 
current is thereafter multiplied by the flicker component 
for voltage in such a way that a product is produced. The 
product is processed in such a way that a flicker power 
n is obtained with a sign value that indicates in which 

25 direction a flickering source is located in relation to 
the measurement point. 

According to an embodiment of the invention, the sign 
value of the flicker power at the measurement point is 
30 negative when the flickering source is located below the 
measurement point and positive when the flickering source 
is located above the measurement point. 




method for deciding 



4 



The method has two important advantages : 




1. The flicker power at the measurement point is 

determined. By this^ means — the flickering source can be 



traced . 



2 . 



The method is also accurate when the measurement 



is carried out via voltage and current transformers 
10 installed in the network. The frequency spectrum of a 
signal with flicker consists of a carrier wave (for 
example 50 or 60 Hz) and of side-bands on both sides of 
the carrier wave, with frequency distance to the carrier 
wave corresponding to the flicker frequencies. The 
15 frequency spectrum of the whole signal package (carrier 
wave and the low-frequency signals) is thus in a narrow 
frequency band around the carrier wave, which means that 
measurement via existing current and voltage transformers 
can be used as these are designed for highest accuracy 
20 around the network frequency. 

| The invention can be carried out by means of a number of 
methods described below. 

25 In a first method (method 1) according to an embodiment 
of the invention, the method is — characterized by includes 
the steps of : 

- recording of an amplitude-modulated current signal i (n) 
comprising signals that originate from the network 
30 frequency f c and the low-frequency amplitude variations 
in the current signal i (n) ; 
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- recording of an amplitude-modulated voltage signal u(n) 
comprising signals that originate from the network 
frequency f c and the low-frequency amplitude variations 
in the voltage signal u(n); 

5 - signal processing of the current signal i (n) in such a 
way that only the low- frequency amplitude variations 
remain in the form of a flicker component for the current 
signal i (n) ; 

- signal processing of the voltage signal u(n) such a way 
10 that only the low-frequency amplitude variations remain 

in the form of a flicker component for the voltage signal 
u(n) ; 

- creation of a product by multiplication of the flicker 
component for current and the flicker component for 

15 voltage; 

- processing of the product in such a way that a flicker 
power 11 is obtained with a sign value that indicates in 
which direction the flickering source is located in 
relation to the measurement point. 

20 

An advantage of the first method is that it does not 
require extensive computer capacity, but can easily be 
implemented in a suitable device. 

25 According to another embodiment of method 1, the method 
comprises : 

the signal processing of the current signal i(n) 
comprising the steps: 

creation of a first demodulated signal by 
30 demodulation of the current signal i (n) ; 

j - filtering off (eliminating) of the signals that 

originate from the network frequency f c in the first 
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demodulated signal in such a way that only the low- 
frequency variations remain in the form of the 
flicker component for current; 
the signal processing of the voltage signal u(n) 
5 comprising the steps: 

creation of a second demodulated signal by 
demodulation of the voltage signal; 
J - filtering off (eliminating) of the signals that 

originate from the network frequency f c in the 
10 second demodulated signal in such a way that only 

the low-frequency variations remain in the form of 
the flicker component for voltage. 

According to another embodiment of method 1, the method 
15 comprises the steps: 

- filtering off of the signals that originate from the 
network frequency f c in the first demodulated signal in 
such a way that only the low-frequency variations 
relating to the flicker component for current remain in 

20 the form of a flicker signal I L F<n) for current; 

- filtering off of the signals that originate from the 
network frequency in the second demodulated signal in 
such a way that the low-frequency variations relating to 
the flicker component for voltage remain in the form of a 

25 flicker signal U LF(n) for voltage; 

- the product creating an instantaneous power signal 
n(n) by multiplication of the flicker signal I L F(n) for 
current and the flicker signal U LF(n) for voltage; 

- the product being processed by the creation of the 
30 average value of the instantaneous power signal ll(n) 

whereby the flicker power n is created with the sign 
value . 
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According to yet another embodiment of method 1, the 
method comprises: 

- the first demodulated signal being created by square 
5 demodulation of the current signal; 

- the second demodulated signal being created by square 
demodulation of the voltage signal. 

According to an embodiment of method 1, the filtering is 
10 carried out with a bandpass filter with a lower limit of 
0.1 Hz and an upper limit of 35 Hz. A preferred upper 
limit is, however, 25 Hz. 



In a second method (method 2) according to an embodiment 
15 of the invention, the method 4?s — characterized by includes 
the steps of : 

- recording of an amplitude-modulated current signal i (n) 
comprising signals that originate from the network 
frequency f c and the low-frequency amplitude variations 

20 in the current signal i (n) ; 

- recording of an amplitude-modulated voltage signal u(n) 
comprising signals that originate from the network 
frequency f c and the low-frequency amplitude variations 
in the voltage signal u(n); 

25 - frequency analysis of the wave form of the voltage 
| signal u(n) by means of an N-point DFT analysis (Discrete 
Fourier Transform) , which gives rise to a voltage vector 
U which contains the frequency spectrum for the voltage 
signal u(n) in the form of N complex voltages; 
30 - frequency analysis of the wave form of the current 
j signal i (n) by means of an N-point DFT analysis (Discrete 
Fourier Transform) , which gives rise to a current vector 
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I which contains the frequency spectrum for the current 
signal i(n) in the form of N complex currents; 

| - creation of a power vector P by means e-£ the 

multiplication, element by element, of the voltage vector 
5 U and the current vector I; 

- multiplication of the power vector P by a weighting 
vector W that eliminates the power component that 
originates from the network frequency, with the power 
vector P comprising partial powers P k concerning power 

10 components from the flickering source, 

- creation of a flicker power n with a sign value by 
| means of summation of the partial powers P k , and 

analysis of the sign value, with the sign value 
indicating in which direction from the measurement point 
15 the flickering source is to be found. 

In the second method, the voltage signal u(n) and the 
| current signal i (n) are signal processed by means — of the 
frequency analysis described above. The product given by 

20 method 1 corresponds to the creation of the power vector 
P in method 2. The flicker components for current and 
voltage in method 1 do not have anything directly 
corresponding to them in method 2, but flicker components 
arise in the form of power components P k concerning power 

25 components from the flickering source after the power 
vector P has been multiplied, element by element, by the 
weighting vector W which comprises the elements w k . The 
processing described in method 1 corresponds to the 
summation of the partial powers P k in method 2 . 

30 

According to an embodiment of method 2, the flicker power 
n is created by means — e-f— the following step: 
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I - summation of the partial powers P k by moans — of- the 
formula : 

An advantage of this method is that no demodulation 
5 residues arise. 

According to another embodiment of method 2, the flicker 
| power FT is created by means of the following steps 

- square demodulation x 2 of the voltage signal u(n); 
10 - square demodulation x 2 of the current signal i (n) ; 

- calculation of the frequency spectrum of the square- 
j demodulated voltage signal by means — an N-point DFT 

analysis (Discrete Fourier Transform) which gives rise to 
the voltage vector (U) ; 
15 - calculation of the frequency spectrum of the square- 
| demodulated current signal by means — e# an N-point DFT 
analysis (Discrete Fourier Transform) which gives rise to 
the current vector (I); 
| - creation of the flicker power n by moans — e# summation 
20 of the partial powers P k which contribute to the flicker 
| phenomenon by moans of the formula: 

where the elements wl k and w2 k replace W k and eliminate 
the power component that originates from the network 
25 frequency and weight the correct amplitudes of the 
frequency components U k and I k/ in accordance with: 
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1 



for \<k<i 



0 fork> i 



w2 k =< 



for \<k<i 



0 for k> i 



where it is assumed that the low-frequency flickers are 
5 to be found in a frequency band up to and including tone 

i (0 < f flicker < i) • 

As shown above, both method 1 and method 2 give rise to 
the flicker power n according to the invention with a 
10 sign value that indicates in which direction a flickering 
source is located in relation to a measurement point. It 
is thus possible to signal process the current signal and 
voltage signal in both the time plane and the frequency 
plane in order to obtain the required flicker power IT. 

15 

The invention also relates to an arrangement comprising 
for carrying out the methods described above. 

According to an embodiment of the invention, the 
20 arrangement comprises means for deciding the direction to 
a flickering source in relation to a measurement point in 
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an electricity network with alternating current with a 
network frequency f c with low-frequency amplitude 
variations from the flickering source. The arrangement is- 
charactcrizod in that it also comprises: 
5 - means for recording an amplitude-modulated current 
signal i (n) comprising signals that originate from the 
network frequency f c and the low-frequency amplitude 
variations in the current signal i(n); 

means for recording an amplitude-modulated voltage 
10 signal u{n) comprising signals that originate from the 
network frequency f c and the low-frequency amplitude 
variations in the voltage signal u(n); 

- means for signal processing the current signal i (n) in 
such a way that only the low-frequency amplitude 

15 variations remain in the form of a flicker component for 
the current signal i (n) ; 

- means for signal processing the voltage signal u(n) in 
such a way that only the low-frequency amplitude 
variations remain in the form of a flicker component for 

20 . the voltage signal u(n); 

- means for creating a product by multiplication of the 
flicker component for current and the flicker component 
for voltage; 

- means for processing the product in such a way that a 
25 flicker power n is obtained with a sign value that 

indicates in which direction the flickering source is 
located in relation to the measurement point. 

According to another embodiment of the invention: 
30 - the means for signal processing of the current signal 
(i (n) comprises : 
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- means for creating a first demodulated signal by 
| mcano of demodulation of the current signal i (n) ; 

- means for filtering off the signals that originate 
from the network frequency f c in the first 

5 demodulated signal in such a way that only the low- 

frequency variations remain in the form of the 
flicker component for current; 

- the means for signal processing of the voltage signal 
u(n) comprises: 

10 - means for creating a second demodulated signal by 

j moans of demodulation of the voltage signal; 

- means for filtering off the signals that originate 
from the network frequency in the second demodulated 
signal in such a way that only the low-frequency 

15 variations remain in the form of the flicker 

component for voltage. 

According to another embodiment of the invention, the 
arrangement comprises: 

20 - means for frequency analysis of the wave form of the 
| voltage signal u(n) by means of an N-point DFT analysis 
(Discrete Fourier Transform), which gives rise to a 
voltage vector U which contains the frequency spectrum 
for the voltage signal u(n) in the form of N complex 

25 voltages; 

- means for frequency analysis of the wave form of the 
| current signal i(n) by mcano — eS an N-point DFT analysis 

{Discrete Fourier Transform) , which gives rise to a 
current vector I which contains the frequency spectrum 
30 for the current signal i(n) in the form of N complex 
currents ; 
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- means for the creation of a power vector P by the 
multiplication, element by element, of the voltage vector 
U and the current vector I; 

- means for the multiplication of the power vector P by a 
5 weighting vector W that eliminates the power component 

that originates from the network frequency, with the 
power vector P comprising partial powers P k concerning 
power components from the flickering source, 

- means for the creation of a flicker power n with a 
10 sign value, by summation of the partial powers P k , and 

- means for analysis of the sign value, with the sign 
value indicating in which direction from the measurement 
point the flickering source is to be found. 

15 BRIEF DESCRIPTION OF DRAWINGS 

The invention will be described below in the form of a 
number of embodiments with reference to . a number of 
figures, in which: 

20 Figure 1 shows an equivalent two-terminal network for an 
electricity network according to the invention; 

Figure 2a shows variations in the effective value U L of 
the voltage; 

25 

Figure 2b shows variations in the effective value I L of 
the current; 

Figure 3 shows a frequency spectrum for an amplitude- 
30 modulated voltage signal with only one low-frequency 
component ; 
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Figure 4 shows a signal flow chart for a measurement 
method according to an embodiment of the invention; 

Figure 5 shows a signal flow chart for a measurement 
5 method according to a second embodiment of the invention; 

Figure 6 shows a signal flow chart for a measurement 
method according to yet another embodiment of the 
invention; 

10 

Figure 7 shows an amplitude characteristic for a bandpass 
filter that is used in the embodiment described in 
association with Figure 6; 

15 Figure 8 shows an outline block diagram of an arrangement 
that can be used for the methods described in the 
different embodiments; 

Figure 9 shows schematically a network comprising a 
20 flickering source Fl, a load LI, a generator G for 
generating an alternating voltage, and 

Figure 10 shows schematically a diagram of the flicker 
power n for a number of sampling points n. 

25 

PREFERRED EMBODIMENTS 

Both the general theory and a number of embodiments will 
be described below. The general theory is necessary for 
an understanding of the embodiments described below. In 
30 the equations that are described, a point between two 
letters in an equation indicates a multiplication of 
vectors, element by element. 
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As mentioned above, the invention relates to a method for 
| deciding — determining in which direction a flickering 
source is located in relation to a measurement point. To 
5 determine the direction of low-frequency variations such 
as flicker, the wave form of the voltage and current must 
be recorded in the phase or phases where the direction is 
to be determined. Thereafter, the recorded information is 
signal processed according to any one of the embodiments 
10 that are described in the following signal flow charts. 
The result of the signal processing is a flicker power II 
with a sign value. The sign value indicates in which 
direction the flickering source is located in relation to 
a measurement point. 

15 

Figure 1 shows an equivalent two-terminal network for an 
electricity network according to previously known 
technology. The electricity network can be divided 
schematically into three parts which are commonly called 

20 generator 1, transmission line 2 and load 3. When a load 
3 is connected, a current I will flow in the circuit. 
This leads to a drop in voltage U Z t arising across the 
internal impedance Z T so that the voltage U L across the 
load drops. If the load is connected in and disconnected 

25 periodically (cyclically) , the current I will also 
increase and decrease cyclically, so that the voltage U L 
decreases and increases cyclically (provided that the 
generator voltage U G remains constant) . 

30 Figure 1 shows a measurement point 17 and a point 18 that 
marks a point above, that is upstream of, the measurement 
point 17, and a point 19 that marks a point below, that 
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I is downstream of, the measurement point 17 — . The terms 
above and below the measurement point 17 are important, 
as it is necessary to indicate in which direction the 
flickering source is located in relation to the 
5 measurement point. 

According to the invention, a recording is taken J^at the 
measurement pointy of an amplitude-modulated current 
signal i (n) comprising signals that originate from a 

10 network frequency f c and low-frequency amplitude 
variations in the current signal i (n) . The low-frequency 
amplitude variations originate from the flickering 
source. In addition, a recording is carried out of an 
amplitude-modulated voltage signal u(n) comprising 

15 signals that originate from a network frequency f c and 
the low-frequency amplitude variations in the voltage 
signal u(n) . Also here, the low-frequency amplitude 
variations originate from the flickering source. 

20 Figures 2a and 2b show outline diagrams of variations in 
the effective value U RM s of the voltage for the voltage U L 
and effective value Irm S of the current for the current 
voltage I. According to the invention, the changes in the 
effective values for current and voltage U RMS and Irms 

25 relate to the instantaneous changes in the voltage and 
| current that arises due to the low-frequency amplitude 
variations from the flickering source. Figures 2a and 2b 
do not show signals that originate from the network 
frequency. It is, however, known that the low-frequency 

30 amplitude variations modulate the network frequency f c . 
Any changes that take place in the current and voltage 
and that originate from the network frequency are 
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negligible. The effective value changes U RMS and I RMS thus 
reflect the instantaneous changes in current and voltage 
that can be derived from the low-frequency variations. In 
Figures 2a and 2b, from the point of view of the 
5 invention, it is thus only of interest to show the 
variations in the effective value Urms of the voltage and 
the effective value I RM s of the current. The variations in 
the effective value U RM s of the voltage and the effective 
value Irms of the current depend, according to the method 
10 according to the invention described below, on whether 
the flickering source is located above or below the 
measurement point 17. 

Figure 2a shows the situation when the flicker power is 
15 propagated from the load towards the generator, with the 
changes in the effective values for current I RM s and 
voltage Urms taking place instantaneously and in antiphase 
(180 degrees out of phase) . When the load increases, the 
current I increases, whereby the drop in voltage across 
20 U Z t increases, which means that U L reduces 
instantaneously . 

Figure 2b shows the situation when the flicker power 
propagates from the generator towards the load, with the 

25 changes in the effective values for current Irm S and 
voltage Urm S taking place instantaneously and in phase. An 
increase of U G gives an instantaneous increase of the 
current I, which gives an instantaneous increase of U L , 
which thus gives a simultaneous change of current and 

30 voltage . 
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In a theoretical consideration of flicker, it is 
expedient to let the current signal and voltage signal 
u{t), i(t) be described as an amplitude modulation. The 
mathematical expression for such signals is given by [1] 
5 and [2], in accordance with: 



The signals consist partly of a carrier wave U c cos (oo c t) 
and I c cos(co c t) respectively, and the low-frequency 
flicker signals U m kcos{co k t) and I mk cos(co k t) respectively. 
10 The index c relates to the carrier wave's component and 
the index m relates to the low- frequency variation's 
component. The index k relates to the index for a 
summation across values k=l to N. 

15 Figure 3 shows a frequency spectrum for an amplitude- 
modulated voltage signal with only one low-frequency 
component (single-tone modulation). The low-frequency 
component relates here to a low-frequency tone that gives 
rise to the variations in amplitude of both current and 

20 voltage. By low-frequency tone is meant here a low- 
frequency signal. 

Figure 3 shows the frequency bands to which the low- 
frequency variation gives rise. In Figure 3, f c relates 
25 to the carrier frequency and f m to the frequency of the 
low-frequency variation. The frequency components with 




[1] 




[2] 
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single-tone modulation are to be found at the frequencies 
f c , f c -f m and f c +fm- 

Figure 3 also shows the amplitude U c for the voltage of 
5 the carrier frequency f c and the amplitude U m /2 for the 
voltage component from the low-frequency variation f c . 

Figure 3 shows that a frequency spectrum of the 
modulating tone forms an upper side-band f c +f m and a lower 

10 side-band f c -f m with half the original amplitude U ra in 
each side-band and located at a frequency distance from 
the carrier frequency f c corresponding to the modulating 
frequency f m . Low frequencies such as the low-frequency 
flicker give rise to side-bands close to the carrier 

15 frequency. The higher the frequency of the flicker, the 
larger the frequency distance of the modulating frequency 
from the carrier frequency. 

The flicker power, IT, is the power that has its origin 
20 in the modulating tones and can be expressed as: 

n = y J^S^- cos K^ + P k ^(jj 7 ™* cos to' + a « >]* = f orth °9° na "y } = 

= Z^-ooWk -*k) = t^*- cosfo ) 
*=i k=\ * 

The formula [3] shows that, after multiplication and 
integration, the individual low- frequency tones in 
25 current and voltage create the flicker power. In order to 
determine this power, the amplitude and phase of the low- 
frequency signals in voltage and current must be known 
and the low-frequency signals must be able to be 
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extracted from the signal package [1] and [2]. This can 
be carried out in several different ways that lead to the 
same result. Examples of different methods are described 
below as different embodiments of the invention. 

5 

Two embodiments will be described below that are based on 
the frequency spectrum of the sampled wave forms u[n] and 
i[n] being determined by carrying out an N-point DFT 
analysis (DFT = Discrete Fourier Transform) . In the 

10 practical case, an FFT analysis ( FFT = Fast Fourier 
Transform) can be used, which is a calculation algorithm 
that provides the same information as the Fourier 
Transform. The output data from the analysis is two 
complex-value voltage and current vectors U and I, which 

15 contain frequency spectrums for u[n] and i [n] in the form 
of k = N complex voltages U k and currents I k . The 
frequency resolution that is obtained depends upon the 
selected sampling frequency (f s ) and the number of 
samples N that are included in the calculation, in 

20 accordance with the relationship Af = f s /N . 

Mathematically, the vectors U and I are expressed: 

= , , c/ 3 3 . . ,] = 0^7, |^ A ,|c/ 2 1 ^ 2 , |t/ 3 1^/? 3 , . . . , |c/^ |^>^ ,] [4] 

/ = [/ lJ / 25 /3,..,^,]=[/ 1 |Za 1 ,|/ 2 |Za 2 ,|7 3 |Z« 3 ,...,|/ w |Za A , 5 ] [5] 

The complex power, S, is created on the basis of the 

25 vectors U and I and contains the power quantities P 

(average power) and Q (reactive power) per tone in 
accordance with: 
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[6] 



The active power Pi is obtained as the real part of Si, 
in accordance with: 



5 The vectors U, I and S contain the frequency spectrum 
from the angular frequency f = 0 Hz to f = f a Hz . There 
is often a desire to weight the information and/or zero 
certain frequencies. A simple and effective way of doing 
this is to introduce the weighting vector W, in 
10 accordance with: 



15 The weighing factor W contains elements, w k , that contain 
constants that are multiplied, element by element, by U, 
I and S. If the voltage vector U is to be zeroed for 
frequencies above half the sampling frequency, the 
following operation is carried out (the point notation in 

20 [9] relates to multiplication element by element) : 



/> = Refo}=Re{i> + jQ,} 



w =[i 



[8] 
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"mod = U ™ 

[9] 



where 



I for \<k< — 

2 I 

0 for K < k<m 



5 U mo d is the modified voltage vector containing the 
frequency spectrum up to half the sampling frequency. 
Above this, the elements of the voltage vector are zero. 

The weighting vector's elements, w k , can be selected so 
10 that the required result is obtained. For example, a 
filter characteristic can be achieved by selecting 
suitable values of the elements w k . In addition, 
different weighting vectors can be introduced for U, I 
and S with the aim of achieving the required result. 



The flicker power, IT, is obtained by summation of the 
powers, P k , which contribute to the flicker phenomenon. 
This is to say: 
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20 



[10] 
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Figure 4 shows a signal flow chart for a measurement 
method according to an embodiment of the invention for a 
first measurement method. The calculation of the flicker 
power n is carried out by frequency analysing the wave 
5 form of the voltage and current. 

The input signals for the signal flow chart in Figure 4 
consist of the two input signal vectors, u[n] and i[n], 
which contain the sampled wave forms for voltage and 
10 current. The distance in time between two indices (n and 
n+1) in the input signal vectors corresponds to l/f 3 , 
where f 3 relates to the sampling frequency. 

By utilizing the previously described theory for the 

15 signal flow chart that is described in Figure 4, the 
flicker power can be calculated. Figure 4 shows two 
blocks designated DFT that refer to a transformation of 
the current signal and voltage signal, i(n) and u(n) 
| respectively, by moans — &f- a so-called Discrete Fourier 

20 Transform. The signal that comes out from the DFT for the 
voltage signal u(n) is designated U and comprises the 
components U k Zp\ which relate to a number of values for 
different angles p k , where the index k indicates which 
point is referred to in the DFT. The signal that comes 

25 out from the DFT for the current signal i (n) is 
designated I and comprises the components IkZcp k which 
relate to a number of values for different angles (p k , 
where the index k indicates which point is referred to in 
the DFT. Figure 4 shows a multiplier 3 in which the 

30 vectors U, I and W are multiplied. The output signal from 
the multiplier is designated P and relates to a power 
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vector comprising components partial powers P k that 
describe the active power in accordance with [7] as the 
real part of the complex power S in accordance with [6] 
and [7] . 

P k = Re{0.5-w k - | U k I Zp k - 1 l k | Z-cpk} 

5 

' [10a] 



Figure 4 shows, in addition, a summation point Z that 
10 indicates the creation of the flicker power II by 
summation of the active powers Pk in accordance with 
[10a] , where : 



15 



n=X Re{0.5-w k -|Uklzp k -|l k |z-(p k } 



[10b] 



The frequency spectrums of the wave forms u[n] and i[n] 
are calculated by an N-point DFT analysis. Thereafter, 
the flicker power n is calculated by utilizing [10] . The 

20 content of the weighting vector W is decisive for the 
final result. Elements in W are to be selected so that 
the power in the carrier wave is zeroed and so that only 
the tones that originate from flicker are included. This 
can be carried out in three ways . The information is used 

25 from either 1) the upper side-band, 2) the lower side- 
band or 3) there is a summation per frequency of the 
information in both side-bands, and the power is 
thereafter calculated by utilizing [10]. 
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In order to exemplify how the calculations are carried 
out, we assume that the carrier wave signal is found in 
element M of the vectors U and I . We also assume that the 
5 low-frequency flicker tones in the vectors U and I are 
found in the elements k with the index: 



10 The number of flicker tones is given by the constant i. 
If the choice is made to utilize the information in the 
lower side-band for the calculation of the flicker power, 
the elements w k in the weighting vector W are selected, 
in accordance with 



The flicker power is thereafter created by utilizing the 
formula [10] . 

20 If the choice is made to utilize the information in the 
upper side-band for the calculation of the flicker power, 
the elements w k in the weighting vector W are to be 
selected, in accordance with 



M-i<k<M and M<k<M + i 



15 
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The flicker power is thereafter created by utilizing the 
formula [10] . 

If the choice is made to utilize the information in both 
side-bands for the calculation of the flicker power, the 
5 elements w k in the weighting vector W are selected, in 
accordance with 




fO for \<k<M-i 

1 for M-i<k<M and M<k<M + i 
0 for k = M 



j The flicker power EI is thereafter created by moans — eS 
the formula 

11 = X Re |^ W k ■ ■ ( U 'M-i+k) + U iM+i-k) )■ • (J r M . i+k) + / (M+ i-k) Y j 

10 

[10c] 

Figure 5 shows a signal flow chart for a second 
15 measurement method according to another embodiment of the 
invention, where the calculation of the flicker power is 
| carried out by means — ef- square demodulation of voltage 
and current, after which a DFT analysis of the low- 
frequency signals is carried out. 

20 

The method is reminiscent of the first measurement method 
described in association with Figure 4, with the 
difference that, before the frequency spectrum is 
calculated, the signals are square demodulated, which is 
25 designated X 2 in Figure 5. By this means, the low- 
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frequency signals are separated from the carrier wave, in 
accordance with the expression below: 



u 2 (0 = \u c + f^U ml . cos(<V + /?, ) j cos 2 (co c t + &) = 

= |f// + 2 • t/ c • cos(o t r +A) + ^j[X cos^/ + p k )J j • I (l + cos(2<y c / + 2fi c )) = 

= ■ ^-+^-2-C/ c g<7„ It cos(<y,/+/? t ) + i •^gC/ mt cos(a> t / + /? fr )j •+ 
+ j^y- + ^ 2 U c- i U »* cos (<V + A ) + j ■ ^2^-* C0S (<V + A )j | • ^ cos(2<y f / + 2yf? e ) 

[11] 

Square demodulation means that two signal packages are 
5 created, separated in frequency. One signal package 
consists of a direct component, the modulating 
frequencies and the modulating signals squared. The 
latter is an unwanted combination product. The second 
signal package contains the same terms but centred in 

10 frequency around double the carrier wave frequency. For 
calculation of the flicker power, only the terms are 
included that are marked by double underlining in [11] . 
One term contains the low-frequency flicker tones 
multiplied by the carrier wave amplitude and the other 

15 term is the modulated signal squared. In terms of size, 
the first term is much larger than the second, which 
means that the squared term affects the result very 
little. Corresponding expressions are obtained when 
squaring the current signal i(t). 

20 

The input signals for the signal flow chart in Figure 5 
consist of the two input signal vectors, u[n] and i [n] , 
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which contain the sampled wave forms for voltage and 
current. The distance in time between two indices (n and 
n+1) in the input signal vectors corresponds to l/f s . In 
X 2 each element in the input signal vectors is squared 
5 and constitutes input data for the N-point Discrete 
Fourier Transform, DFT. The output signal from the 
respective DFT is complex-value vectors U and I 
containing frequency spectrums for u[n] and i[n] (both 
amplitude information and phase information) with a 
10 frequency resolution of Af = f s /N. The content in U and I 
is expressed by: 



\^..J N ]=\l x \Za x \l 2 \Za 2 ,\l,\Za,,...,\l N \Za N ] 



The 



complex power, S, is calculated by 



5 




15 | 



[12] 



The 
the 



sought flicker power n is then obtained by 



formula 




20 
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[13] 



The elements wl k and w2 k zero the frequency components 
that do not cause flicker and weight the correct 
amplitudes of the frequency components U k and I k , in 
accordance with 



U- for \<k<i\ 



0 for k>i 



I 

10 [14] 



f— for \<k<i\ 
0 for k>i \ 



[15] 



15 In [13], [14] and [15], it is assumed that flicker tones 
are found in a frequency band up to and including tone i 
(0 < f fucker ^ i) • The correct value for the weighting 
factors can be identified by studying [11] . The method is 
not limited to the weighting factors that are given 

20 above, other constants in the weighting vector being able 
to be used to obtain the required filter effect. 
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Figure 6 shows a signal flow chart according to yet 
another embodiment of the invention for a third 
measurement method, where the flicker power is calculated 
by square demodulation of voltage u(n) and current i (n) 
5 where the low-frequency tones are filtered by the use of 
bandpass filters 2A and 2B respectively. 

Instead of using DFT analysis and weighting vectors as 
described in the two previous embodiments described in 
10 association with Figures 4 and 5, in this measurement 
method bandpass filters 2A, 2B are used for filtering the 
low-frequency flicker tones. 

The input signal vectors u[n] and i [n] contain the 
15 sampled wave forms for voltage and current. In the 
components 1A and IB, the carrier wave (signal from the 
network frequency) and the low-frequency flicker tones 
are separated, similarly to the measurement method 
described in association with Figure 5, by each sample 
20 being squared. Only the low-frequency flicker tones are 
allowed to pass through the bandpass filters 2A and 2B. 
As input signals to the multiplier 4, there are thus only 
the low-frequency flicker tones in voltage u LF [n] and 
current i LF [n] . The output signal from the multiplier 4 
25 is the instantaneous flicker power ri[n] = u LF [n] x i LF [n]. 
The flicker power II is obtained by the integrator 5 
creating an average value for the instantaneous power 
EI(n). This can be carried out with a digital filter in 
the form of, for example, a low-pass filter. 

30 

The bandpass filters 2A and 2B are dimensioned to obtain 
a lower limit frequency of 0.1 Hz and an upper limit 
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frequency of 25 Hz. Alternatively, the bandpass filters 
can be used that are defined in the standard IEC 61000-4- 
15 that describes a flicker algorithm. 

5 The transmission function for the filter in the IEC 
standard is expressed as 



10 



H(s) = ^Q) x s \+sla 2 

s 2 +2As + e» 2 (s + l/ry 3 )0 + l/<y 4 ) 



[16] 



The coefficients in [16] are to have values in accordance 
with the following table. 



k = 1.74802 


X = 271-4. 05981 


co 1 =27T-9. 154 94 


co 3 =27t-l .22 535 


(02=271-2.2797 9 


co 4 =27t-21 . 9 



15 Table 1. 



Figure 7 shows the amplitude characteristic for the 
bandpass filter with transmission function in accordance 
with [16] and coefficients in accordance with Table 1. 

It is also possible to select bandpass filters with a 
different characteristic to that described in [16]. For 
example, an Mth order Butterworth or Chebishev filter can 
be selected. 
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Figure 8 shows a block diagram of the hardware for the 
instrument. The instrument is constructed around a signal 
processor 7 which administers the measurement information 
and carries out necessary calculations according to the 
5 measurement methods that are described in the Figures 4-6 
and 9. The signal processor 7 also controls the sampling 
process for the conversion of analogue signals to digital 
signals (the A/D conversion) . The recorded signals-^ that 
is_j_ the wave form of current and voltage, are obtained 

10 either from current transformers and voltage transformers 
in the network or from measurement sensors that are to be 
found in the instrument. The signal processor can be a 
computer or a logic circuit, or some other suitable 
device for controlling devices and for signal processing 

15 of signals. 

Figure 8 shows a signal-conditioning device 8 for the 
recorded voltage signal. The measured voltage is signal- 

| matched by m e ans of a resistive voltage division, whereby 

20 the correct input signal level is obtained for the 
following step which is an anti-alias filter 10. 

Figure 8 also shows a signal-conditioning device 9 for 
the recorded current signal. The signal level for the 

25 current channel is matched to the instrument, either via 
a low-ohm shunt, with the voltage drop across the shunt 
being amplified and becoming the input signal level for 
the following step which is an anti-alias filter 11, or 
alternatively, the current signals can be obtained from 

30 current tongs connected to the instrument. 
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• The task of the anti-alias filters 10, 11 is to prevent 
weighting distortion which arises if the recorded signal 
has a frequency content that exceeds half the sampling 
frequency (see theory for the sampling theorem). The 
5 anti-alias filters can be implemented in the form of an 
analogue Sallen-Key low-pass filter according to Figure 7 
and have a limit frequency corresponding to half the 
sampling frequency. 

10 The level-matched and filtered signals are sampled in 
sampling devices 12, 13 with a sampling frequency, for 
example 6400 Hz. The digital raw data in the form of 
sampled amplitude values is saved in a measurement memory 
14 and later constitutes the input data for the 

15 measurement methods described above. 

Software that controls the signal processor 7 is found in 
a program memory 15. The finished result, that is the 
flicker power with sign value and internal impedance, can 
20 be shown both numerically and graphically in a display 
device 16. The display device can be any known device for 
numerical and graphical display, for example a VDU . 

Figure 9 shows schematically a network comprising a 
25 flickering source Fl, a load LI, a generator G for 
generating an alternating voltage. In Figure 9, the 
current direction for the current I is shown by an arrow 
in the connecting lines. The current goes from the 
generator G to the flickering source Fl and the load LI. 
30 Figure 9 also shows a first measurement point Ml. In 
association with Ml, there is marked a point above 18 the 
measurement point Ml and a point below 19 the measurement 
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point Ml. Figure 9 also shows a second measurement point 
M2 . In association with M2 there is marked a point above 
18 the measurement point M2 and a point below 19 the 
measurement point M2 . The flickering source Fl emits into 
5 the network a low-frequency amplitude variation that 
propagates in the direction shown by a solid-line arrow 
20 in Figure 9. 

In the method according to the invention described above, 
10 the flicker power n indicates in which direction the 
flickering source is located in relation to a measurement 
point. According to an embodiment of the invention, the 
sign value of the flicker power is negative when the 
flickering source is located below 19 the measurement 
15 point and- positive when the flickering source is located 
above 18 the measurement point. 

At the first measurement point Ml, a negative sign is 
thus obtained for the flicker power n, as the flickering 
20 source Fl is located below 19 the first measurement 
point. This is because the low- frequency variations in 
current and voltage are in antiphase at the first 
measurement point Ml . 

25 At the second measurement point M2, on the other hand, a 
positive sign is obtained for the flicker power n, as 
the flickering source Fl is located above 18 the second 
measurement point M2 . This is because the low-frequency 
variations in current and voltage are in phase at the 

30 second measurement point M2 . 
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Figure 10 shows schematically a diagram of the flicker 
power n for a number of sampling points n during a 
certain period of time when recordings were made of the 
modulated current signals and voltage signals i (n) , u(n). 

5 

Figure 10 shows a first curve Kl on the negative lower 
part of the diagram. The diagram also shows a second 
curve K2 on the upper part of the diagram. The first 
curve Kl corresponds to a power signal which, after the 

10 creation of the average value, gives rise to a flicker 
power n with a negative value and thus corresponds to 
the flicker power n that was obtained at the first 
measurement point Ml in Figure 9. The second curve K2 
corresponds to a power signal which, after the creation 

15 of the average value, gives rise to a flicker power n 
with a positive value and thus corresponds to the flicker 
power n which was obtained at the second measurement 
point M2 in Figure 9. 

20 Figure 10 also shows that Kl inverted about the sampling 
axis n corresponds to K2 . This has been demonstrated in 
experiments and the explanation is that the flicker power 
n changes sign when the propagation of the low-frequency 
flicker is changed from going against the power direction 

25 of the basic tone to going with the power direction of 
the basic tone, or vice versa. It is thus necessary to 
know in which direction the generator and load 
respectively are to be found in relation to the 
measurement point in order for it to be possible to 

30 interpret the sign value in the correct way. The part of 
the diagram that shows the period between zero and where 
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Kl and K2 respectively start, shows a period of time when 
the flickering source Fl is not connected in. 
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